Unconventional superfluid order in the F-band of a bipartite optical square lattice 



Matthias Olschlager, Georg Wirth, and Andreas Hcmmcrich 
Institut fur Laser-Physik, Universitat Hamburg, Luruper Chaussee 149, 22761 Hamburg, Germany 

(Dated: August 30, 2010) 

We report on the first observation of bosons condensed into the energy minima of an F-band of a 
bipartite square optical lattice. Momentum spectra indicate that a truly complex-valued staggered 
angular momentum superfluid order is established. The corresponding wave function is composed 



of alternating local F 2x 



iF, 



2y3-3y 



■orbits and local S'-orbits residing in the deep and shallow 



wells of the lattice, which are arranged as the black and white areas of a checkerboard. A pattern 
of staggered vortical currents arises, which breaks time reversal symmetry and the translational 
symmetry of the lattice potential. We have measured the populations of higher order Bragg peaks 
in the momentum spectra for varying relative depths of the shallow and deep lattice wells and find 
remarkable agreement with band calculations. 

PACS numbers: 03.75.Lm, 03.75.Hh, 03.75.Nt 



The physics of condensed matter systems is often de- 
termined by electrons, which possess orbital degrees of 
freedom involving an intrinsic anisotropy due to multiple 
orbital orientations. Orbital physics plays a central role 
for magnetism, superconductivity, and transport prop- 
erties of rare earth and transition metal compounds in- 
cluding high T c cuprate systems or heavy fermion sys- 
tems, which are subject of intense research since more 
than two decades (JHI- The advent of optical lattices 
(i.e. quantum gases arranged in synthetic lattices formed 
by light) has raised hopes that certain aspects of such 
systems could be studied in a precisely controlled envi- 
ronment without many of the complexities usually asso- 
ciated with material systems 0, la]. Unfortunately, the 
wave function of bosons in their ground state according 
to a consideration by Feynman, referred to as the "no 
node theorem", is positive definite under very general 
circumstances [g, 0] • This appears to significantly limit 
the use of bosons for simulating many-body systems of 
interest. Fermions, on the other hand, are significantly 
harder to prepare in optical lattices [a Q and the real- 
ization of large filling factors required to access relevant 
orbital physics appears difficult, if not impossible. It has 
long been speculated, that the use of higher bands could 
be a way to explore the interplay between superfluidity 

The devel- 



and orbital physics also with bosons 
opment of clever new techniques has made it possible to 
populate higher bands I n Ref- 12 1 atoms could 



be excited into the P-band of a two-dimensional (2D) 
optical lattice and short-lived (< 1 ms) coherence was 
established along specific directions. Cross-dimensional 
coherence could be established more recently in a bipar- 
tite lattice due to cross-dimensional tunneling junctions 
realized by local S'-orbits in every second well [19( . How- 
ever, the .D-and F-bands in optical lattices have remained 
practically unexplored by experiments. 

In this article we report the first observation of bosons 
condensed in the P-band of a quasi 2D bipartite optical 
lattice. Full cross-dimensional coherence with a life-time 



on the order of 10 ms is established. The observed mo- 
mentum spectra exhibit a characteristic pattern of sharp 
maxima, which are well explained by a complex-valued 
superfluid order parameter composed of alternating local 
F2x s -3x + * p2y 3 -3y" orrj its and local S'-orbits. A pattern 
of staggered local angular momenta and staggered vorti- 
cal currents arises commensurable with the plaquettes of 
the lattice with the consequence of broken time-reversal 
symmetry. The proposed nature of the superfluid order is 
confirmed by evaluating the Bloch functions correspond- 
ing to the observed condensation quasi-momenta, which 
are derived via two-dimensional band calculations. We 
measured the population ratios between higher order and 
lowest order Bragg peaks in the momentum spectrum for 
varying relative depths of the shallow and deep lattice 
wells and find remarkable agreement with calculations 
based upon these Bloch functions. 

By crossing two optical standing waves derived from 
laser beams with 100 ^m l/e 2 -radius and a wavelength 
A = 1064 nm we produce a (quasi 2D) light shift potential 
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providing two classes of (tube-shaped) lattice sites (de- 
noted as A and B) arranged as shown in Fig. 1(a). Here, 
k = 2tt/X, rj is experimentally adjustable (around unity) 
and e is fixed to sa 0.9 due to imperfect reflection optics 
used in the experiment. For details we refer to Ref . (l9| . 
Adjustment of the parameter 8 lets us tune the difference 
of the well depths between A and £>-sites. For 6 < ir/2 
the .4-sites are more shallow than the S-sites and vice 
versa. A change of 9 by ir/2 can be obtained in less than 
0.2 ms. 

Efficient population of excited bands is obtained by 
optimizing a population swapping procedure described 
in more detail in Ref.[l3|. Initially a Bose-Einstein con- 
densate (BEC) of Rubidium ( 87 Rb) atoms is prepared 
and the lattice potential is ramped up within 80 ms to 
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FIG. 1: (a) The bipartite lattice comprises two classes of 
lattice sites denoted by A and B. The grey area shows 
the Wigner Seitz unit cell of the „4-sublattice. (b), (c) and 
(d) Observed populations of Brillouin zones (BZs) after ap- 
plying the population swapping procedure with final values 
of (9 /it, V /E icc ) = (0.61, 5.8), (0.66, 8.3), (0.69, 13.3), respec- 
tively, (e) Schematic of the 6th (blue and purple areas) and 
7th (orange areas) BZs. The eight red circles mark the points, 
where the atoms are observed to gather in (c). In the cen- 
ter of the plot the first BZ is reconstructed by translations of 
subsets of the 6th BZ via reciprocal lattice vectors indicated 
by the red dashed arrows, (f) 3D representation of the image 
in (c). (g) energy surface of the 7th band with degenerate 
minima indicated by black disks. 



Vo/E Tec = 16.6 for a value 9 < tt/2 such that the well 
depth of the B-sites significantly exceeds that of the A- 
sites (E lec = h 2 k 2 /2m denotes the recoil energy with the 
atomic mass m). A ground state lattice is thus formed 
with most atoms residing in the deeper ZJ-wclls. The 
large well depth yields nearly complete suppression of 
tunneling. Subsequently 9 is rapidly changed (within 0.2 
ms, which is shorter than the nearest neighbor tunneling 
time) to a hnal value Of above 7r/2 such that now the A- 
wells are significantly deeper than the £>-wells. Finally, 
Vo is adiabatically decreased during 0.6 ms to admit tun- 
neling again. Optionally a 2 ms long phase is appended, 
where 9 is adiabatically tuned to some desired value fol- 
lowed by a variable hold time. 

In Fig.l(b),(c), and (d) the population of different 
bands is illustrated. To obtain these pictures, popula- 
tion swapping is carried out with 9f and Vo optimized 
for maximal population of the desired band. The lat- 
tice beam intensity is then exponentially decreased with 
a time constant of 430 [is and after 30 ms absorption 
imaging is applied. This yields images of momentum 
space, where the population P[n, q] in the n-th band for 
some quasi-momentum q and energy E[n,q] is mapped 
to some point within the n-th Brillouin zone (BZ) re- 
lated to q by a reciprocal lattice vector. This mapping 
requires that E[n,q] docs not cross some other E[n',q] 
during the adiabatic switch off process. As a conse- 



quence of such band crossings a part of the band pop- 
ulation may be mapped to an adjacent BZ. The image in 
Fig.l(b) (with 9/tt = 0.61 and V /E Tec = 5.8) matches 
well with the 4th BZ, while in Fig.l(d) (with 9/tt = 0.69 
and Vo/E lec = 13.3) the 9th BZ is found to be the most 
populated. This directly indicates the numbers of the 
bands that have been populated to be n = 4 and n = 9, 
respectively. For the case of interest in this article in 
Fig. 1(c), where 9/tt = 0.66 and V /E Tec = 8.3, the as- 
signment of a single BZ is not possible due to a band 
crossing in the adiabatic mapping process between the 
6th and 7th band. Instead the atoms share the 6th and 
7th BZ with emphasis on the 6th BZ. This is illustrated 
by Fig. 1(e), where the shapes of the 6th (blue and pur- 
ple) and 7th (orange) BZs are sketched. One recognizes 
in Fig.l(b),(c), and (d) that the BZs are not evenly pop- 
ulated but rather a large fraction of the atoms resides at 
well localized momenta, which for the case of Fig. 1(c) are 
highlighted by the eight red circles at the inner boundary 
of the 6th BZ in Fig. 1(e). This localization in momen- 
tum space is more easily observed in a 3D representa- 
tion, which for the case of Fig. 1(c) is shown in Fig. 1(f). 
As is substantiated below, a collision driven condensa- 
tion process sets in on a timescale comparable to that 
of the band population procedure, which acts to redis- 
tribute the atoms into the energy minima of their band. 
Collision aided band decay, repopulating the 1st BZ, is 
only observed on a significantly longer timescale. 

In the remainder of this article we focus on the 
situation in Fig. 1(c), which corresponds to the 7th 
band (the lowest of the four possible F-bands). In 
Fig. 1(g) we show a plot of the 7th band derived from 
a band calculation involving a Fourier expansion of the 
Bloch- functions with 11 harmonics in each dimension 
and the potential of Eq.(l). The parameters 9/tt = 
0.66 and Vo/E rcc = 8.3 correspond to the final set- 
tings in Fig. 1(c). Indicated in Fig. 1(g) by black disks, 
four local minima of the energy surface arise at the 
quasi-momenta K(x,i), K(_x,-i), K(i i), K(_i i x), where 
K( I/;J[1 ) = i/ifc [vx + iMy) with integers v, \i and x, y denot- 
ing the unit vectors in x- and y-directions. Condensation 
is expected to occur at these points, which are found to 
coincide in their energies to better than 10 -3 -E roc for a 
wide range of settings of e and r\ in Eq.(l). In particular, 
a local imbalance of the available standing wave intensi- 
ties, which is unavoidable due to the use of finite sized 
Gaussian beams, does not notably lift the degeneracy. 
Our band calculation shows that, when Vq is ramped 
to zero, a crossing between the 7th and the 6th band 
occurs for quasi-momenta in the vicinity of the conden- 
sation points K( ±1 ,±1). Thus, our experimental band 
mapping procedure used to obtain the picture in Fig. 1(c) 
maps these points into the 6th rather than into the 7th 
BZ. This is shown in the center of the BZ illustration in 
Fig. 1(e) where the 6th BZ is mapped onto the first BZ 
via translations with reciprocal lattice vectors (indicated 
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FIG. 2: In (a), (b), and (c) momentum spectra are shown for 
a hold time of 1 ms, Vo/E Iec = 8.3 and 0/tt = 0.66, 0.70, 0.75, 
respectively. Observations (calculations) are shown in the up- 
per (lower) row. (d) Population ratio for higher order and 
zero order Bragg peaks plotted versus 9. Observations (cal- 
culations) are shown by the symbols (solid lines), (e) Tem- 
poral evolution of the populations of Bragg peaks. The solid 
lines are exponential fits applied in the wings of the graphs to 
determine relaxations times, (f) Temporal evolution of band 
populations for Vo/E r<sc — 8.3 and 9/tt — 0.67. The solid lines 
are exponential fits, (g) 1/e-times for depopulation of the 6th 
and 7th BZ and repopulation of the 1st BZ plotted versus 8. 



by the red dashed arrows). Note, that the points of in- 
creased population in the measured BZ-plot in Fig. 1(c) 
(corresponding to the red circles in Fig. 1(e)) are mapped 
into the condensation points of Fig. 1(g) (black disks). 

In Fig. 2(a), (b) and (c) momentum spectra are shown 
for a hold time in the lattice of 1 ms, a well depth parame- 
ter V /E rcc = 8.3 and 0/tt = 0.66, 0.70, 0.75, respectively. 
The upper row shows the experimental observations ob- 
tained by rapidly (< 1/xs) disabling the lattice potential 
and absorption imaging after 30 ms of ballistic expan- 
sion. The presence of sharp Bragg peaks at quasi mo- 
menta K(„ M ) for odd integers v, fj, clearly demonstrates 
cross-dimensional coherence. One recognizes the absence 
of a zero momentum component (y = fi = 0). In fact, the 



lowest order components arise at the four quasi-momcnta 
K(±i.±i) that we already identified as the condensation 
points in Fig. 1(c). In the lower row the observations are 
contrasted with calculations of the Fourier components of 
the coherent superposition Hfp = 0k (1 1} +* </>k (1 _ 13 of the 
real-valued Bloch functions corresponding to the two in- 
equivalent condensation points. The specific choice of the 
relative phase "i" will be justified below. A remarkable 
agreement with the observations arises with regard to the 
8 dependence, which is studied more quantitatively in 
Fig. 2(d). In this graph the symbols show the population 
ratios between higher order ((i/, fi) — (3, 1), (3, 3), (7, 1)) 
Bragg peaks and the zero order Bragg peak ((v,fJ,) = 
(1,1)) observed for varying values of 9. The peak pop- 
ulations are obtained by counting the number of atoms 
in a small circular region covering an individual Bragg 
peak and subtracting the number of atoms found in a 
surrounding ring-shaped area of the same size. The the- 
oretical curves (solid lines) do not involve any free pa- 
rameters. In Fig. 2(e) we study the timescales for the 
formation and decay of coherence by evaluating momen- 
tum spectra at fixed values Vo/E lec = 8.3 and 9/tt = 0.66 
for different holding times. The graph shows the Bragg 
peak populations (obtained as in (d)) for the peaks iden- 
tified in the inset. During the first few ms, which re- 
flects the collisional timescale for the condensation pro- 
cess and the tunneling time, all peak populations increase 
before decay sets in with a timescale corresponding to 
the collisional relaxation of the band populations. Both 
timescales are determined by exponential fits as 1.8 ms 
and 9.1 ms for the zero order peak (black diamonds) and 
2.9 ms and 10.1 ms for the first order peak identified by 
red rectangles. 

We also investigated the temporal evolution of the 
band populations. To this end the 7th band was pop- 
ulated as explained above and after a variable hold time 
band mapping was applied in order to produce images 
as in Fig. 1(c). The populations found in the differ- 
ent BZs are plotted versus the hold time in Fig. 2(f) for 
V /E rcc = 8.3 and 9/tt = 0.67. The 4th and the 5th BZ, 
which maintain nearly constant population during the 
observation time, have been omitted for better legibility. 
The graph shows a complex population redistribution, at- 
tributed to collisional dynamics. As a consequence, the 
initially prepared population in the 6th and 7th BZ de- 
cays and the 1st BZ is repopulatcd. The corresponding 
relaxation times are determined by means of exponential 
fits (solid lines) to be 12.3 ms and 19.9 ms, respectively. 
Similarly, the relaxation times for other values of 9 are 
determined and plotted in Fig. 2(g). The difference of the 
relaxation times reflects the fact that the 1st BZ is mostly 
populated via transfer from lower lying bands with life 
times longer than that of the 7th band. One recognizes, 
that around 9/tt sa 0.68 the relaxation times attain max- 
imal values. The local S'-orbits of the wave function 'Sp 
(a detailed discussion of its geometry follows below)) be- 
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FIG. 3: Orbit configuration of order parameter ^f. The grey 
areas characterize the antinodal structure of the orbits, the 
colored numbers indicate the local phases. 



maximally isotropic shape, such that the atoms (which 
interact repulsively) can best avoid each other. We thus 
conclude that, in presence of repulsive collisions, ^J/f rep- 
resents the true ground state of the F-band. For P-orbits 
a similar prediction is discussed in Refs. [13, ED]. The 

as- 
sumption of an incoherent mixture of phases ^>k (1 ±) and 
^k (1 _!) would necessarily imply spatial separation, since 
at the shallow wells both Bloch functions share equivalent 
local S-orbits. A phase separation scenario, however, re- 
quires excess kinetic energy at the phase boundaries due 
to inhibited tunneling and excess mean field energy. A 
compensating local anisotropy, sufficient to enforce local 
condensation in a single condensation point, is not avail- 
able in our experiment. Phase separation thus appears 
energetically unfavorable. 



come maximally developed for this setting of 9. Because 
the ground state wave function has vanishing amplitude 
in the shallow wells, its overlap with "Fp is reduced for 
this case, and thus collisional decay into the ground state 
is inhibited. 

The wave function 'Fp employed in the calculations 
of the momentum spectra has remarkable properties il- 
lustrated in Fig. 3, where its local amplitude (grey areas 
indicate large amplitude) and phase (colored numbers) 
is sketched. In the deep wells 'Ff closely resembles the 
superposition V'p.o] i i ^[0,3] °f eigenfunctions ^[n.m] of 
a 2D harmonic oscillator with n, m oscillator quanta in 
x- and y-directions. The spatial shape of ^f around 
the center of the deep wells is thus approximated by 
(2x 3 — 3x)+i (2y 3 — 3y). In contrast, in the shallow wells, 
"Fp mimics the harmonic oscillator S-orbit ipm^]. The 
checkerboard-like arrangement of S-orbits and F-orbits 
with alternating angular momentum provides equal lo- 
cal phases on both sides of the tunneling junctions and 
thus maximizes the tunneling efficiency. The inner and 
outer dashed rectangles denote the unit cells of the lat- 
tice and of ^f, respectively. The sublattice of shallow 
wells possesses a pattern of staggered vortical currents 
commensurate with its plaquette structure, which match 
with the alternating orbital currents in the deep wells. 
As a consequence <Ff breaks the translation symmetry 
of the lattice and time-reversal symmetry. Order param- 
eters with similar properties have been recently predicted 
in the ground states of driven optical lattices [2C 
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tions <^k (1>1) and 
thus a wave function </>k 



The observed condensation momenta imply that the 
underlying state should be approximated by a coherent 
superposition or an incoherent mixture of the Bloch func- 
k m and 0K a _!)■ If coherence is assumed, and 
j r c0k (1 _ 1 . with some com- 
plex constant c, immediately c ~ i follows in order to 
reproduce the observed momentum spectra. The setting 
c = i maximizes the local angular momentum in the deep 
wells and thus minimizes the mean field energy. The 
reason is that in this case the local F-orbits acquire a 
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